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Abstract. Using Wilson's numerical renormalization group (NRG) technique 
we compute zero-bias conductance and various correlation functions of a double 
quantum dot (DQD) system. We present different regimes within a phase diagram 
of the DQD system. By introducing a negative Hubbard U on one of the quantum 
dots, we simulate the effect of electron-phonon coupling and explore the properties 
of the coexisting spin and charge Kondo state. In a triple quantum dot (TQD) 
system a multi-stage Kondo effect appears where localized moments on quantum 
dots are screened successively at exponentially distinct Kondo temperatures. 



PACS numbers: 72.10.Fk, 72.15.Qm, 73.63.Kv 

1. Introduction 

The development of nanotcchnology has stimulated studies of transport through 
coupled quantum dots systems, where at very low temperatures Kondo physics as well 
as magnetic interactions play an important role. A double-dot system represents the 
simplest generalization of a single-dot system which has been extensively studied in the 
past. Recent experiments demonstrate that similar realistic devices can be constructed 
HIE.! 4], which enables direct experimental investigations of the competition between 
the Kondo effect and the exchange interaction between localized moments on the dots. 
One manifestation of this competition is a two stage Kondo effect that has recently 
been predicted in multilevel quantum dot systems with explicit exchange interaction 
coupled to one or two conduction channels [5, 6J. Experimentally, it manifests itself 
as a sharp drop in the conductance vs. gate voltage G(Vg) [Z] or as non-monotonic 
dependence of the differential conductance vs. drain-source voltage dI/dVd s (Vd s ) [5]. 

Another line of research of correlation effects in quantum dots concerns the 
influence of localized phonon modes on electron transport through quantum dots 
and molecules. The influence of a strong, localized electron-phonon interaction [5] 
can be simulated by the introduction of an effective negative-C/ Anderson model. In 
such system the charge Kondo effect leads to screening of charge fluctuations on the 
localized impurity by low- frequency pair fluctuations in the leads |10j . It has recently 
been shown that in the case of large electron-phonon coupling the original single- 
impurity Anderson model, coupled to local phonon degrees of freedom, can be mapped 
onto an anisotropic Kondo model [llj . Increasing the electron-phonon coupling leads 
to a suppression of the conductance plateau width [TTJ [U] as a function of the gate 
voltage. 
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We study a double quantum dot (DQD) in a side-coupled configuration (figure[T]), 
connected to a single conduction-electron channel. Similar systems have been studied 
using various analytical techniques such as the non-crossing approximation [12], 
embedding technique [13j or slave-boson mean field theory |14[ I15j , and recently also 
using more accurate numerical renormalization group (NRG) calculations [16] . 

In this work we study three distinct cases. In the case, already studied in detail 
in Refs [El \17\ [18] , when effective Coulomb interactions on DQD are positive and 
the intra-dot overlap is large, we find wide regimes of enhanced conductance as a 
function of gate-voltage at low temperatures due to the Kondo effect. Regimes of 
enhanced conductance are separated by regimes where localized spins on DQD are 
antifcrromagnetically (AFM) coupled. In the limit when the dot a is only weakly 
coupled, the system enters the "two stage" Kondo regime [El US], where we again 
find a wide regime of enhanced conductivity under the condition that the high- and 
the low- Kondo temperatures (Tk and T K respectively) are well separated and the 
temperature of the system T is in the interval T K <C T <C Tk- 

In the second case we explore the coexistence of the spin and charge Kondo 
state with different Kondo temperatures that are realized when there is an effective 
attractive interaction on one of the quantum dots. We show that spin and charge 
susceptibilities are screened at the corresponding spin and charge Kondo temperatures. 
Such a spin-charge Kondo state can be realized when two molecules, one with a strong 
phonon mode and the other with a localized orbital where electrons experience strong 
Coulomb interaction, would be coupled to metallic leads. 

Finally, in the case of three side-coupled quantum dots we demonstrate that 
a multi-stage Kondo effect is realized provided that exchange couplings between 
successive quantum dots are smaller than the corresponding Kondo temperatures. 




Figure 1. (Color online) Side-coupled configuration of quantum dots 



2. Model and method 

The Hamiltonian that we study reads 

H = S d (n d - 1) + S a {n a - 1) - (4<V + a U*) (!) 

(7 

+ ^(n d -lf + ^(n a -lf (2) 

+ Y ekC l* Ck ° + Y Vd< ^> (4<A + d Uk<r) > (3) 
ha ha 

where n d — ^ CT d\d a and n a = a a a a- Operators d\ and a\ are creation operators 
for an electron with spin a on site d or a. On-site energies of the dots arc defined by 
e = 5 — {7/2. For simplicity, we choose the on-site energies and Coulomb interactions 
to be equal on both dots, S a = S d — S. Coupling between the dots is described by 
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the inter-dot tunneling coupling t^. Dot d couples to both leads with equal hopping 
t'. As it couples only to symmetric combinations of the states from the left and the 
right lead, we have used a unitary transformation [2U] to describe the system as a 
variety of the single-channel, two-impurity Anderson model. Operator c\ a creates a 
conduction band electron with momentum k, spin a and energy efe = —D cos k, where 
D = 2t is the half-bandwidth. The momentum-dependent hybridization function is 
Vd(k) — — (2/\/N + T)t' sin k, where N in the normalization factor is the number of 
conduction band states. 

We use Meir-Wingreen's formula for conductance in the case of proportionate 
coupling [21j which is known to apply under very general conditions (for example, the 
system need not be in a Fermi-liquid ground state) with spectral functions obtained 
using the NRG technique [22} [23] [24] [25] . At zero temperature, the conductance is 

G = G irTA d (0), (4) 

where Go = 2e 2 /h, Ad(u>) is the local density of states of electrons on site d and 
T/D = {t'/tf. 

We use the density-matrix [55] [T7] version of the standard NRG method [52] . In 
all our NRG calculations we took into account the rotational invariance in the spin 
space as well as in the isospin space, i.e. SU{2) sp - l - a ® SU (2)i sosvm . 



3. Phase diagram of the DQD system 

As a part of the introduction, we first review some findings on the model with repulsive 
Hubbard interaction Ud — U a = U > as already given in Refs. [17l[T8]. We first 
present the phase diagram of the DQD system as depicted in figure [2] Regime of high 
conductance is represented by the dark-gray area. High conductance is a consequence 
of the Kondo regime where (n) ~ 1. Full lines defining the Kondo regime are given by 
the following analytical expressions: Si = td(2^/l + (U/4td) 2 — 1) and 62 = (U/2 + td). 
A spin-singlet regime forms in the strong coupling regime as soon as (n) ~ 2, while in 
the weak coupling regime, i.e. for tg/D < 1 the spin-singlet state persists down to the 
following condition J e gp = Tk- In the opposite case, when J ff < Tk, exponentially 
small scale, the two-stage Kondo regime appears for T < T^. At low temperatures, 
the conductance is enhanced when < T < J c g. 

We now explore in detail the limit when td 0. In this case one naively expects 
to obtain essentially identical results as in the single-dot case. It turns out, however, 
that a non-zero inter-dot coupling, td ^ 0, represents a relevant perturbation to the 
system in the limit T — > 0. A clear physical picture is obtained by comparing the two 
relevant energy scales: the exchange interaction J c g and the single impurity Kondo 
temperature. In the region defined by J e g < Tk and S < U/2, the system enters the 
two-stage Kondo regime. In this regime, a gap of width T K opens in the density of 
states Ad(u>) [16] which consequently leads to a decrease of the conductance when the 
temperature is below T K . The lower Kondo temperature T K depends exponentially 
on T K as [HI [19] 

T° ~T K eM~uT K /J cS ). (5) 

When by decreasing 5 or by increasing td, the exchange interaction overcomes Tk, i.e. 
J e ff > Tk, the gap in Ad(u>) is given by J c ff. In this regime the two local moments on 
DQD form a spin-singlet. 
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Figure 2. Phase diagram of a DQD for U/D = 1, obtained using analytical 
estimates as given in the text. Grey area represents the region of the Kondo 
regime where S ~ 1/2, (n) ~ 1, and G/Go ~ 1. In the shaded area, called 
spin - singlet regime, where S ~ and (n) ~ 2, spin-spin correlation function is 
enhanced, i.e. (S a • S<j) < —0.5 and G/Go ~ 0. The two-stage Kondo regime is 
explained further in the text. 



The two-stage Kondo regime can be clearly observed in Fig [3^ where we show 
the impurity part of the spin susceptibility Xs 

Xs(T) = ^ «S Z 2 ) - (S 2 Z ) ) . (6) 

The first expectation value in this expression refers to the system with the double 
quantum dot, while the second refers to the system without the dots. The two- 
stage Kondo effect manifests as two successive decreases of the susceptibility, first at 
T ~ Tjc ~ 1CT 7 followed by the second at T ~ Tg- [16]. In this case, the two-stage 
Kondo effect occurs for ta/D < 1.6 10~ 4 . 

The same effect can also be seen in the temperature-dependent impurity 
contribution to the entropy S(T), defined via 

sm.2zfi.2i2i, p) 

where E = (H) = Tr (He~ H ^ kBT ^>) and F = -k B T hxTr (e- H / kBT ). At T - T K the 
entropy first decreases from 2 In 2 to In 2, see Fig[3jD. This drop indicates the Kondo 
screening of the local moment on the dot-<i while the moment on the dot-a remains 
unscreened and nearly free. The second drop takes place at T ~ T^. 

With increasing 5, Tk increases, and the condition for the two-stage Kondo effect 
is satisfied at increasingly larger td/D. In Fig 2] we present Xs(T) and S(T) calculated 
at 5/D — 0.4. The two-stage Kondo effect sets in at substantially larger value oitd/D 
than in the 6 — case, i.e. at tg/D ~ 0.01. 

In Fig [S] we explore how the two Kondo temperatures reflect in the shape of the 
frequency dependent spectral function Ad(ui), computed on the dot d. At u ~ U/2 all 
curves display a charge-transfer peak. At lower frequencies spectral functions again 
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Figure 3. Temperature-dependent spin-susceptibility a), and entropy b) for 
various values of td/D, computed at 5 = 0, Y/D = 0.03, U/D = 1. See also 
Ref. [16]. 




Figure 4. The same as in Fig|3] except for 8 = 0.4. 
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increase around u> ~ Tx , and then decrease around w ~ T&-. All calculations were 
done at the effective temperature T/D ~ 10 -16 , which is lower than T^- for all t^, 
except for t^/D ~ 10~ 8 . Taking into account the Meir-Wingreen formula [3T] for the 
zero-bias conductance, equation (g]), we realise that in the limit T — > 0, G/Go — > for 
any finite i d /L>. At small but finite T, such that Tg- < T < !>, the gap in = 0) 

closes and the conductance approaches its unitary limit G/Go ~ 1. 




to/D 



Figure 5. Normalized spectral function A^(lo) of DQD system. System 
parameters are identical to those used in Fig [3] See also Ref. 1 161 ■ 



4. Spin versus charge Kondo effect 

In this section we investigate a system with Ud = — U a = U > 0. The negative U on 
dot-a simulates the effect of strong electron-phonon interaction. Beside the impurity 
spin susceptibility and entropy we also compute the normalized charge susceptibility, 
defined as 

Xc(T) = ^«/ z 2 )-(/ 2 2 )o), (8) 

where I z is the z-component of the total isospin of the system. The isospin operators 
are given by 

*i = ^2vl a <raa>Vi,a', (9) 

ota' 

where the Nambu spinor nj on the impurity orbitals are defined by 




Isospin raising operator is defined as / + = djdj — ajaj. Similarly, isospin operators can 
be defined on the Wilson's chain and the total isospin operator can then be defined by 
a sum of isospin operators over all orbitals of the problem (impurities and conduction 
band). For S — 0, both I 2 and I z commute with H and / and I z are additional good 
quantum numbers. 

Since the system consists of two quantum dots, one with positive- and the other 
with negative- U we expect that the dot-d would display a spin Kondo effect where 
the local spin moment would be screened by spin fluctuation in the lead, while dot-a 
would display a charge Kondo effect where charge degrees of freedom on dot-a would 
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be screened by pair fluctuations in the leads. Due to different overlaps with the leads 
we expect the spin Kondo temperature Tks to be different from the charge Kondo 
temperature, Tkc- 

In figure [5] we present results of Xs{T), Xc(T) and S(T) for various inter-dot 
coupling strengths. At large td/U > 0.3 both susceptibilities as well as the entropy 
show a rapid drop with the temperature, characteristic for a system where the ground 
state of the DQD is a spin and isospin singlet ( state with no spin nor charge 
fluctuations). This finding in agreement with the solution of the DQD system in 
the atomic limit. For td > \^3U/6 ~ 0.29C/, the ground state is nondegenerate spin 
and isospin singlet (S = 0, 1 = 0) with energy 

E (S = 0, 1 = 0) = -2t d . (11) 

In the opposite case, the ground state becomes four-fold degenerate spin doublet 
(S = 1/2), isospin doublet (/ = 1/2) with energy 

E (S = 1/2, 1 = 1/2) = - ^jt\ + E/ 2 /4. (12) 

In the td ^ limit, the latter state is composed of a spin (up or down) moment 
on the dot-d and zero (I z = —1/2) or doubly occupied (I z = +1/2) dot-a. All 
curves in figure [5] for td/D < 0.22 show two separate Kondo screenings: spin Kondo 
effect is visible in Xs{T) at Tks and charge Kondo effect is visible in Xc{T) at 
Tkc < Tks- This is corroborated by the shape of S(T) that shows a plateau at 
2 In 2 at T > Tks, m agreement with the four-fold degenerated 5 = 1=1/2 ground 
state, a drop at T ~ Tks, followed by a plateau at In 2 due to the uncompensated 
isospin moment (degenerate zero and doubly occupied states on the dot-a) , and finally 
a drop at T ~ Tkc- Decreasing td/D affects Tkc much more than Tks- A simple 
explanation where the effective hybridization of the dot-a strongly depends on td/U 
while hybridization of the dot-d becomes i^-independent for small td seems to be 
adequate. 

In Fig. [7| we present spectral functions of the DQD system. In all figures there 
is a charge-transfer peak, located at oj = y/t d + U 2 /4. There is an additional high- 
frequency peak which is, for td < \/3U/6, located at u = + U 2 /4 — 2td- In 
this regime, i.e. for td < V3U/6, a sharp increase in Ad{ui) at ui <~ Tks signals 
development of the spin Kondo peak on the dot-d, while at lower frequency uj ~ Tkc 
a gap opens in Ad(u>). In contrast, there is an additional increase of spectral weight in 
A a (uj) below this frequency. In the limit of lo — > there is an isospin Kondo peak in 
A a (uj), mirrored by a gap of width Alu ~ 22V c within the spin Kondo peak in Ad(ui). 

5. Three-stage Kondo model in a side-coupled TQD model 

A natural question arises when studying the two-stage Kondo effect in a side-coupled 
quantum dot: could one realise a multi-stage Kondo effect in a series of side-coupled 
quantum dots. We consider the following model of three impurities, one coupled 
directly to a continuum and two coupled in chain to the first one. Instead of solving 
the more realistic three Anderson impurities model as depicted in figure [5] we solve a 
simplified version of three Kondo impurities: 

H = -ffband + 2_j Hi + H c + H\2Z (13) 
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Figure 6. Temperature-dependent XsiXc and S for various values of t^/D, 
computed at S = 0, T/D = 0.03, U d /D = -U a /D = 1. 




Figure 7. Spectral functions A^Ilo) a) and c) A a (uj) b) and d) of DQD system. 
System parameters are identical to those used in Fig [6] 
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Figure 8. Side-coupled configuration of three quantum dots. 



#band = ^ e kc\ a Cko, (14) 
ho 

Hi = S(n i -l) + ^-(n i -lf, (15) 



2 

-j=EM c Ldi*+^c*«r), (16) 

kcr 

H 12 3 = JiSi ■ S 2 + J 2 S 2 • S 3> (17) 

where Ji = — 4if 2 /^ an d -^2 = — 4t| 3/^- In. the case when J\ < Tk, where Tk is 
the single-impurity Kondo temperature of the dot-1, and J 2 < T K where T K is the 
second-stage Kondo temperature, yet another low-energy scale 

T K ~ exp(-aT° / J 2 ) (18) 

appears in the problem. The system under the above conditions undergoes three 
S=l/2 Kondo screening cross-overs, i.e. a three-stage Kondo effect, as can be seen in 
Fig.! 

At T > Tk, there is the local moment regime where S/kg ~ 3 In 2 and the local 
moment approaches Tximp/(g^B) 2 ~ 3/4, characteristic of three noninteracting local 
spins. For temperatures below the first Kondo temperature, i.e. T < Tk, the local 
moment on dot-1 couples into a singlet with conducting electrons while spins on dots-2 
and 3 remain uncoupled, which gives S/Ub ~ 2 In 2 and X Ximp / (gf-B) 2 ~ 1/2. There 
is a second drop of S and Ximp at T < T K where the spin on the dot- 2 becomes 
screened by the quasiparticles in the leads forming a Fcrmi-liquid with the screened 
moment on the dot-1. The leftover of the entropy S/Ub ~ In 2 as well as the local 
moment T 'xi mp / (g^B) 2 ~ 1/2 in this temperature range are due to unscreened spin 
on the dot-3. Finally, there is a total screening of the local moment at T < T K . 

The numerical renormalization group eigenvalue flow is shown in Fig. 1101 Each 
Kondo screening induces an additional 7r/2 quasiparticle scattering phase shift while 
one "free" local moment is simultaneously removed from the problem. Each 7r/2 
phase shift corresponds to a change of the boundary condition for the conduction 
band electrons from periodic to anti-periodic (or vice-versa) or, from another point 
of view, each Kondo screening removes one conduction band site from the Wilson 
chain |22j . As a consequence, we observe an alternation between two sets of fixed 
point energies, each set corresponding to one possible boundary condition. Note that 
neither quantum numbers nor the degeneracy of the levels composing one energy level 
remain the same after two such alternations: this is due to the successive removal of the 
"free" local moments from the problem. The degeneracy of the group of states having 
(approximately) the same energy decreases after each Kondo cross-over; in particular, 
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Figure 9. Temperature-dependent susceptibility and entropy for a side-coupled 
chain of three Kondo-like impurities. 



the eigenvalues corresponding to high total spin are pushed to higher energies. Of 
course, on the final strong-coupling fixed point is stable. 

6. Conclusions 

We have explored different regimes of the side-coupled DQD. In the case of positive 
Hubbard interactions when quantum dots are strongly coupled, wide regions of 
enhanced, nearly unitary conductance exist due to the underlying Kondo physics. 
Results were summarised in a phase diagram of the DQD structure. When DQD is 
doubly occupied, conductance is zero due to formation of the spin-singlet state which 
is effectively decoupled from the leads. When quantum dots are weakly coupled unitary 
conductance exists at finite temperatures when two electrons occupy DQD due to a 
two-stage Kondo effect as long as the temperature of the system is well below Tk and 
above T%, 

It should be noted that if the signs of both electron-electron interaction parameters 
Ud and U a are inverted, we obtain an equivalent problem, the only change being that 
the roles of spin and isospin are permuted. If both quantum dots have negative U, 
we will thus observe enhanced conductance due to the charge Kondo effect for strong 
inter-dot tunneling coupling and suitable S, and two-stage charge Kondo effect for 
weak inter-dot coupling. 
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Iteration number N 

Figure 10. Numerical renormalization group eigenvalue flows for odd and even- 
length Wilson chains of the side-coupled triple quantum dot system. The system 
exhibits three-stage Kondo screening at three different Kondo temperatures. The 
levels are labeled by the quantum numbers / and S, the total isospin and total 
spin of the system. 



We have explored a system where one of the dots is strongly coupled to phonon 
degrees of freedom. This effect was simulated by the introduction of an effective 
negative-// interaction on one of the dots, while U on the other dot was kept positive. 
The main goal was to explore the possibility of the coexistence of the spin and charge 
Kondo effect. We found out that, indeed, spin and charge Kondo effect can coexist. 
Even though our model contains two leads, due to mirror symmetry the underlying 
problem falls into the class of the single-channel Kondo problems. It is worth stressing 
that as a consequence of the spin-charge separation in the leads, a single channel is 
sufficient to completely screen two different moments - spin as well as isospin. For 
Ud > 0, U a < and small tj, the local moment (spin) on the directly coupled dot 
occurs at a higher spin Kondo temperature, while the isospin degree of freedom on the 
side-coupled dot is screened at much lower charge Kondo temperature. For Ud < 
and U a > the roles of the spin and charge are simply reversed: there is a charge 
Kondo effect on the directly coupled dot and spin Kondo effect on the side-coupled 
dot. 

In the case of TQD, a multi-stage Kondo effect exists for the appropriate choice 
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of the exchange coupling constants. In the zero temperature limit, the conductance 
approaches unitary limit in the particle-hole symmetric case, since the TQD fixed 
point is identical to the single-impurity fixed point with unitary conductance in the 
particle-hole symmetric point. In fact, this conclusion can be generalized in two ways. 
First of all, the fixed point and the quasiparticle phase shift in the case of three 
impurities is the same irrespective of the values of exchange constants. Depending on 
the exchange constants, we can either have three-stage Kondo effect (weak exchange 
coupling), S=l/2 Kondo effect for an antiferromagnetic chain composed of the three 
impurities (both exchange constants large), or two of the neighboring impurities couple 
into a local singlet, while the spin on the third impurity is Kondo screened (if one of 
the exchange constants is much larger than the other). The conductance is unitary 
in all these cases. We can also generalize our findings for the case of more than 
three impurities. In the case of N side-coupled quantum dots we predict unitary zero 
temperature conductance in the case of odd N and zero conductance for even N. For 
suitable exchange constants, the screening can occur in several (N or less) stages and 
the conductance as a function of the temperature will be non- monotonous. 
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